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Abstract 

The important consequences of the recent results of the numerical evaluations of eighth and tenth 
order QED contributions to the anomalous magnetic moment of electron are commented. The 
correctness of the results of the numerical evaluation of new eighth order QED corrections to the 
electron anomaly are supported by the demonstration of their consistency with the new analytical 
expressions for the QED contributions to a e from the diagrams with fourth-order light-by-light 
scattering muon and tau-lepton loops. The consistency of the similar results are demonstrated in 
the case of eighth order massive dependent contribution to the muon anomalous magnetic moment, 
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One of the most precise at present experimental results in the modern particle physics is 
the measurement of the electron anomaly a e = (g e — 2)/2 [l|, j2[, which gives 

a e = 1159652180.73(0.28) x 10~ 12 [0.2Appb] . (1) 

The previous stage of the Standard Model theoretical prediction 

af = aj ED + a e (hadrons) + a e (weak) (2) 

was summarized in the most detailed review on the subject 0]. The perturbative QED 
contribution to Eq.([2]) is defined as 

af ED = Ai + A 2 {m e /m II ) + A 2 {m e /m T ) + A 3 (m e /m M , m e /m r ) . (3) 

Up to recently theoretical and phenomenological applications were based on the following 
expression for its dominant term Ai 

5 2 

A x = V Af ) = 0.5 (-) - 0.32847896557919378 ...(-) (4) 

l=1 7T \7T/ \7T/ 

+ 1.181241456587... (^) 3 - 1.9144(35) + 0.0(3.8) 5 

where / is the number of loops of the Feynman diagrams, which are contributing to the 
corresponding perturbative QED expression. Three first coefficients, presented in Eq.(jl]) in 
the numerical form, were evaluated analytically. The first term was calculated by Schwinger 
Q, the second correction was evaluated by Petermann jEf and Sommerfield p]. This re- 
sult was confirmed later on by Terentiev [7J, who used different technique. The project of 
analytical evaluation of all three-loop QED contributions to a e was completed by Laporta 
and Remiddi js|] . The cited value for the four- loop QED correction to a e was numerically 
obtained in Refs.jsj, [Io[ by Kinoshita and collaborators. The rough CODATA estimate of 



the coefficient of the 5- loop term in Eq. <^j, namely A^ = ±2\A^\ [11], gave the idea what 
might be theoretical uncertainties in the value of a e due to unknown up to recently total value 
of the tenth-order QED effects. The project of their direct numerical evaluation, started in 
2005 by Kinoshita and Nio in Ref. [l2j, continued in the series of works of Refs. [l3j]- (2l|. 
was successfully completed in May of 2012 by Ayoama, Hayakawa, Kinoshita and Nio [221 ] . 

In general, the coefficient A^ is defined by the contribution of 12672 diagrams, which 
were classified into 32 gauge-invariant subsets. The result of evaluation of the final, most 
complicated Set V, together with more detailed numerical calculation of the eighth-order 



contribution to a e , was reported recently in Ref. 22]. The long-expected expression for the 
tenth-order massless contribution [22| is : 

A[ W) = 9.16(58) (5) 



while at the eight-order level the following new results were obtained [22J: 

Af = -1.9106(20) (6) 
a ( § ) = -1.9097(20) . (7) 

The c4 8 )-term differs from Eq.(ll]) due to the inclusion of complete mass- dependent contri- 
butions, numerically evaluated in the process of the works of Refs.[l3[-j2l[. These massive- 
dependent effects, summarized in in the numerical form in Ref. [22| . read 

A { 2 8 \m e /m^ = 9.222(66) x 10" 4 (8) 
A ( i\mJm T ) = 8.24(12) x 10~ 6 (9) 
Af\m e /m^m e /m T ) = 7.465(18) x 10~ 7 (10) 
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While obtaining these results the CODATA-2010 report [23| mass ratios m e /m M = 
4.83633166(12) x 10" 3 , m e /m T = 2.87592(26) x 10" 4 , mJm T = 5.94649(54) x 10" 2 , with 
the fixed value of r-lepton pole mass m T = 1776.82(16) MeV 1 , were used. The results of 
Eq.()H]), Eq.flH]) and Eq.f llOp are presented in more detailed form in Table I of Ref.j22[ for 12 
gauge-invariant groups of 4-loop diagrams, contributing to A 2 and A 3 -terms of Eq.([n]). The 
summary of massless and massive-dependent contributions to and A ( 2 w \m e /m lJ ) are 



presented in Table II of Ref. 22] . Note, that the sum of 10-th order massive dependent terms 
result in the small value of the overall massive correction A { 2 W \m e /m^) = -0.00382(39) Q. 



It should be stressed, that the determination of the concrete value of A^ 10 ''- and a^-terms 



is of real importance. Indeed, prior the work of Ref. |22| the scientific community faced with 
the unique case, when the comparison of the low-energy experimental result of Eq.([T]) with 
the massless perturbative QED predictions of Eq.(T4]), supplemented with analytically known 
from the works of Refs. jlBf- (3l[ massive-dependent forth and sixth-order contributions into 
A 2 and yl 3 -terms and with the well-known values for a e {hadrons) and a e {weak) (see the 
review j3[), namely 

a e (weak) = 0.0297(5) x 10~ 12 (11) 
a e (hadrons) = 1.671(19) x 10~ 12 , (12) 

led to the value of the inverse fine coupling constant 

a' 1 = 137.03599084(33) (39) (13) 

with theoretical uncertainties ±39 x 10~ 8 , comparable with experimental ones ±33 x 10~ 8 . 

The dominant contribution to the theoretical error, namely ±30 x 10~ 18 , was defined by 
the CODATA estimate of the coefficient of the 0(a 5 )-term in Eq.(|4]). The additional sizable 
uncertainty in Eq.( fl3|) came from the theoretical error of the numerical evaluation of the 
eighth-order contributions to Eq.(J4]) performed in Refs.jsf, [Io[. There were no additional 
theoretical errors in the sum of cited above analytically evaluated forth and sixth-order 
massive- dependent contributions A 2 and A 3 -terms. Their total contribution can be extracted 



from the summary part of Ref. [3l| and reads: 



a\ 2 _„_ R /a x3 



A 2 (m e /m ll ) = 5.19738667(26) x HT 7 - - 7.37394155(27) x 10~ 6 - (14) 



7T/ \7T 

A 2 (m e /m T ) = 1.83798(34) x 10~ 9 (^ 2 - 6.5830(11) x 10~ 8 (^ 3 (15) 

A 3 {m e /m^m e /m T ) = 0.1909(1) x 10~ 12 . (16) 

The errors in Eq. ( }T4"j) -Eq. ffl6|) are related to the indicated above uncertainties of the 
CODATA-2010 values of the ratios of leptons masses. 

As to the uncertainties of the important contributions from a e (hadrons) and a e (weak), 
they were not sensitive at the previous stage of comparing theoretical and experimental 
predictions for a e 



The results of Ref. 22] and Eq. (jSJ allowed to solve this intriguing problem and to make 



theoretical uncertainties in the analog of Eq. ([T3"j) less important , than experimental ones. 



1 Recently measured value of the r-lepton mass m T = 1776. 69^q \\l ± 0.15 MeV [24| should not change a lot 
the uncertainties of these ratios. 
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Indeed, substituting new eighth- and tenth-order QED effects of Ref. [22J into the the 
procedure of the comparison with the precise experimental result of Eq.flT]), the authors of 
Ref. [22j obtained more precise value of the inverse fine coupling constant: 

a' 1 = 137.0359991657(68) (46) (24) (331) [0.25ppb] . (17) 

Here the first and second errors are related to the uncertainties of the numerical evaluation of 
the eighth and tenth-order QED corrections. The third error is determined by the combined 
uncertainties of the hadronic and electroweak contributions to a e , which start to manifest 
themselves at this more precise level of perturbative QED calculations, while the fourth and 
hugest uncertainty is determined by the experimental error in Eq.([T]). 

In view of the importance of the results of the complicated eight and tenth-order numerical 
calculations it is highly desirable to perform the independent cross-checks if not all of them, 
but at least of some their parts. 

In this work this problem is studied by analyzing the consistency of the numerical results 
for the massive -dependent contributions into A < 2\m e /m fl ) and A2 8 ^(m e /m T )-corrections to 
the electron anomaly, which are described by the subset of diagrams, formed by external 
light-by-light scattering muon and r-lepton subgraphs with extra virtual photon, propagat- 
ing inside this subgraph. The corresponding analytical expressions for the leading term of 
heavy-mass expansions of these contributions follow from the obtained in Ref.[32| result of 
analytical QCD calculations of the hadronic light-by- light scattering contribution to the lep- 
ton anomalous magnetic moments a e and with hadronic effects, modeled by the internal 
light-by-light-scattering quark loop, crossed by virtual gluon. Taking into account that in 
QED one has Cp = 1 and ajjs = a(l + 0(a 2 )) (see Ref.(33| for details) we get the cor- 
responding analytical contribution of the leading term in heavy lepton mass expansion of 
the the eighth-order QED correction to a e from this subset of massive- dependent light-by- 
light scattering graphs. The results contain the contributions of the Riemann ^-functions 
Ck = E^i(l/^) k and polylogarithmic functions a k = Li k (l/2) = Eklil 1 / 2 ^) an d read: 

,(8),,, „, Arr ^,/a\ 4 / 473 A , 2 312 > , , 42853 A 5 771 , 

Ar(X 1+i . Ibl, NLO) - = C2 hi 2 2 + (2 hi 3 2 C A + Ca n 2 

2 v + ' ' V 180 s 405 s 2880 s 360 s 

473 , 4 52 , - 8477 473 416 34727, 23567 \ 2 „ /a\ A 

+ In 4 2 In 5 2 + a 4 + a 5 + Cs Cs Y X* - 

1080 675 2700 45 45 2400 s 1440 S / fr[ 1+1 \n J 

where Xi+i = (mi/mi +! ) with i — 1, 2 are the ratios of the electron mass mi = m e and the 
muon mass 777,2 = Trip or r-lepton mass = m T . The abbreviations in the parenthesis labels 
the expressions for the subsets of graphs, which contain the next-to-leading order (NLO) 
approximation of the light-by- light (Ibl) scattering subgraph. 

Substituting the known values for the transcendental and polylogarithmic functions into 
Eq.flXED we get: 



Af(X 1+i ,lbl,NLO)(^y = 1.77831... £ X 
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The CODATA-2010 values and errors of X\ = m e /rrip and X2 = m e /m T are fixing the 
numerical expressions for the terms we are interested in : 

Af\x 2 ,lbl,NLO)(^y = 4.15948(21) x 10~ 5 (^) 4 (19) 

Af(X 3 ,lbl,NLO)(^\ = 1.47082(26) x 10" 7 (20) 
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The results of the numerical calculations for the similar contributions are presented in 
Table 1 of Ref. 1221 and read 



A ( *\x 2 ,lbl,NLO)(^y = 4.105(93) x 1(T 5 (21) 



Af ) {X^lbl,NLO)[^j = 1.431(95) x 1(T 7 (^j . (22) 

Taking into account the related uncertainties one can observe good agreement with the 
numbers of Eq. ([T9~l) and Eq. (|20p . which follow from analytical expression of Eq. ffTS]) . This can 
be considered as the strong check of the results of computer numerical calculations of eight- 
order contributions to a e , summarized in Ref. 22|. Moreover, minor discrepancy between 
Eq.( )20|) and Eq.( )22|) may indicate, that in the case of larger lepton mass (m r ) the exact 
numerical results of Ref.[22[ may be sensitive to taking into account higher terms of large 
mass expansion of the corresponding eight-order light-by-light-type diagrams, contributing 
to a e . 

Let us have a look whether the similar feature is manifesting itself in the case of the 
comparison of the results of analytical and numerical calculation of the similar eight-order 
light-by-light scattering corrections to the anomalous magnetic moment of muon a M . In this 
case the expression , analogous to Eq. lTTBl . reads 

Af \X 4 Jbl,NLO)(^y = 1.77831... X%(^y . (23) 

Using the CODATA-10 value for X4 = m M /m r we get 



A { £\x^lbl,NLO)(^y = 6.2888(11) x 10 _3 (^) 4 • (24) 
This number is in good agreement with the result of the numerical eighth-order computer 



calculations, obtained in Ref. |34| in the process of completing numerical evaluation of 



tenth-order QED contributions to the muon anomaly a^, namely with the result 



A^>{X^lbl,NLO)[-) = 6.106(31) x 10 3 f — V . (25) 



7T/ \7T 



However, as in the case of r-lepton light-by-light-scattering eighth order contributions to 
a e , the numerical comparison of the analytically-based result of Eq. ( 1241) with the result of 
numerical calculations of Eq.f l25p seem to indicate the sensitivity to still unknown higher 
terms of large mass expansions of the analytically evaluated massive-dependent Feynman 
graphs. 

In any case our considerations demonstrates the reliability of the definite results of the 



complicated important numerical QED calculatiobns from Ref. [22| and Ref.[34|. More de- 
tailed considerations of new 10-th QED order results for using the renormalizat ion-group 



inspired studies of Ref.[35[ may be done in future. 

This work grew up from the detailed studies of the results for the eighth- and tenth-order 
QED contributions to the anomalous magnetic moment of electron and muon, started in 
January of 2012 during the visit to Theory Division of CERN. I am grateful to the members 
of Ph-Th Division of CERN for hospitality. It is the pleasure to thank A.E. Dorokhov, M. 
Nio, M.Passera and T. Kinoshita for discussions of various scientific subjects, related to this 
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